Laser-to-proton energy transfer efficiency in laser-plasma 

interactions 



E. Fourkal, I. Velchev, and C-M. Ma 
Department of Radiation Physics, Laser Laboratory, 
Fox Chase Cancer Center, Philadelphia, PA 19111, U.S.A. 
(Dated: June 10, 2008) 

Abstract 

It is shown that the energy of protons accelerated in laser-matter interaction experiments may 
be significantly increased through the process of splitting the incoming laser pulse into multiple 
interaction stages of equal intensity. From a thermodynamic point of view, the splitting procedure 
can be viewed as an effective way of increasing the efficiency of energy transfer from the laser light 
to protons, which peaks for processes having the least amount of entropy gain. It is predicted that 
it should be possible to achieve > 100% increase in the energy efficiency in a six-stage laser proton 
accelerator compared to a single laser-target interaction scheme. 

PACS numbers: 52.38.-r, 52.38.Kd, 52.38.Ph 
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I. INTRODUCTION 



Ion acceleration by high-power lasers has attracted significant attention in recent years 
from the scientific community due to its potential applications in different branches of physics 
and technology. The physical characteristics of accelerated protons such as high coUimation 
and high particle flux make them very attractive for possible applications in controlled 
nuclear fusion 1, 2!, material science and hadron therapy [4, 5]. 

The physical processes responsible for ion acceleration during laser-matter interaction 
are understood on a qualitative level. The initially proposed theoretical model for ion 
acceleration at the back surface of the target heavily relied on the notion of quasi-neutral 
plasma expansion into vacuum . In this model the accelerating electric fleld is generated 
due to space-charge separation in a narrow layer at the front of the expanding plasma cloud, 
which is assumed to be neutral. In the interaction of ultrashort and ultra-intense laser 
pulse with a solid structure, the assumption of quasi-neutrality must be abandoned. The 
results of computer simulations j^, suggest that the interaction of petawatt laser pulses 
with plasma foils leads to the formation of extended regions where plasma quasi-neutrality 
is violated, a factor that has to be taken into account when considering ion acceleration by 
ultra-intense pulses (J ~ 10^^ W/cm^). According to the model, the incoming laser pulse 
quickly ionizes the target pushing some of the electrons out of it through the action of the 
ponderomotive force. A strong electrostatic fleld (~ TV/m) is set up between the expanding 
electrons and the target that fleld ionizes the thin hydrogen-rich layer present at its back 



surface. Subsequently, the protons are accelerated in this electrostatic field (10|, For 
thicker targets (> 2 /im) a shock wave acceleration mechanism has also been proposed 12| 
in which a laser acts as a piston driving a fiow of ions into the target and launching an 
electrostatic shock at the front of the target with high Mach number M=fshock/c— 0.2-0.3. 
Protons, refiected off the shock front may get accelerated to velocities up to fions = 2f shock- 
In a recent multi-parametric particle-in-cell (PIC) simulation study of interaction between 
a clean (no prepulse present) high-power laser pulse and thin double-layer target |13|] the 
authors have mapped maximum proton energy regions as functions of target electron density 
and its thickness as well as laser pulse length for different laser intensities and spot sizes. 
According to the scaling laws obtained by the authors, in order to accelerate protons to the 
energy range of a few hundred MeV (e.g. required for hadron therapy applications where 
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protons with energy 250 MeV can reach any disease site throughout the patient's body), 
one needs to pump a few hundred joules of energy or equivalently several tens of petawatt 
of power (for laser pulse duration Lp ~100 fs) into the laser pulse, provided that optimal 
conditions for the laser/target parameters are met. Currently available lasers, specifically 
compact table-top systems, operating in the sub-picosecond regime provide energy on the 
order of (Si ~10 J. According to the scaling laws obtained in Ref. 13j this energy is insufficient 
to accelerate protons to the required therapeutic energy range of 200-250 MeV. Therefore, it 
would be interesting to investigate whether there is a way to increase the maximum proton 
energy or equivalently the efficiency of energy transfer from the laser pulse into accelerated 
protons without increasing the energy of the laser pulse. 

In a recent work [1^], we have shown that in a double layer target system, the acceleration 
conditions for protons are far from the optimal due to the fact that protons are expelled from 
the back surface of the substrate before the maximum electric field is established and as a 
result experience reduced acceleration potential. It has been shown that by modifying the 
proton dynamics (through splitting the pulse into two interaction stages), higher final proton 
energies can be achieved. Up to 30 % increase in the final proton energy as compared to a 
single interaction stage was predicted through the double splitting procedure. The natural 
extension of the idea presented above would be to find out whether the energy transfer 
efficiency from the laser pulse to protons can be further improved. An obvious question 
that one may ask is: if a single interaction splitting leads to ~ 30 % increase in the proton 
average energy, then how would the introduction of n interaction stages infiuence the final 
proton energy? 

In this work we will show that splitting a single interaction scheme into n stages gradually 
increases the energy transfered from the laser pulse to a proton beam with each additional 
splitting, thus increasing the final energy of the proton beam. We will offer a thermodynamic 
analogy to this effect, which helps elucidate the underlying physics. In fact, the problem 
at hand is related to the question of what the most efficient way of transferring the energy 
from a hot object (laser) to a cold (protons) one is. It is a known fact that any energy 
exchange process between the constituents of a closed system is accompanied by an increase 
in the total entropy of the system, no matter what the interaction is (the second law of 
thermodynamics). The magnitude of the entropy increase however depends on the manner 
in which the energy exchange occurs and smaller increment in the entropy yields higher 
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energy transfer efficiency [15|. According to classical thermodynamics, a process occurring 
near equilibrium is the most efficient in terms of generating the maximum amount of work 
and the least amount of entropy. A reversible thermodynamic process merely needs to be 
slow enough for all thermally-interacting constituents to equilibrate. This means that the 
high efficiency of energy conversion can be attained, provided that the relaxation time is 
much shorter than the characteristic time of the process. We show that the most efficient 
heat exchange occurs when the cold and hot objects are split into n equal pieces and each 
individual hot piece is put into thermal contact with each individual cold piece (without 
mixing them) in a sequential manner (see Appendix). In the end, initially hot/cold pieces 
are put back together to form a new cold/hot object correspondingly. As the number of 
splits increases, the entropy change AS* for the whole process decreases and in the limit 
n oo, AS 16|. In this asymptotic case the initially cold object becomes hot (with 



temperature equal to the initial temperature of the hot object) and initially hot object 
becomes cold (with temperature equal to the initial temperature of the cold object) and the 
perfect (completely reversible) heat exchange process is achieved. It should be noted that 
the introduced heat exchange problem does not violate the second law of thermodynamics, 
since the energy always flows from the hot object to the cold and increased heat transfer 
efficiency is achieved through the reduction of the amount of energy exchanged in each 
individual interaction step, which ultimately reduces the overall energy "wasted" during the 
heat flow. Therefore, it is reasonable to conjecture that when the laser pulse is split into 
n sub-pulses of equal intensity lo/n that are made to interact with n targets (laser-target 
multi-system constitutes a closed system), the energy transfer efficiency (kinetic energy of 
the accelerated protons) increases due to a decrease in the total entropy gain for the closed 
system (laser light, all particles and flelds). Just as in the case with hot/cold reservoirs, the 
splitting procedure is an effective way of reducing the energy wasted during the laser-target 
interaction, thus increasing the energy transferred from the laser pulse to protons (adiabatic 
acceleration). 

Bringing the entropy argument into the description of laser-mater interaction physics 



essentially nonequilibrium process) and giving it a prominent role comes naturally a: 
ong string of successful applications of the minimum entropy production principle 
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ter 



13, 



19| to various systems including plasmas. The qualitative analogy we offer is a further 



demonstration of the far reaching consequences of this principle that may eventually be used 
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in a quantitative analysis of the physics of laser-matter interactions. 



II. MULTI-STAGE PROTON ACCELERATION IN 2D PARTICLE-IN-CELL SIM- 
ULATIONS AND 3D MODEL 

In order to study how the introduction of multiple interaction stages influences the final 



proton energy, we have used 2D PIC simulations [20|, |2l| to model the interaction between 
the laser pulse and several targets. The calculations were performed in a 2048 x 1024 
simulation box with a grid size A = 0.04 fim and total number of simulated quasi-particles 
2.6 X 10^. Periodic boundary conditions for particles and electromagnetic fields have been 
used. The initial conditions were chosen to correspond to realistic experimental parameters, 
where linearly polarized {p polarization) relativistically intense (/o=1.92x 10^^ W/cm^, A= 
800 nm), ultrashort (Lp ^ 30 fs) Gaussian laser pulse (focal spot size D=3.2 fim at FWHM) 
is normally incident in a Cu target of thickness 400 nm and transverse dimension 5 /im. 
The electron density as well as the ion charge state in the target are ne=3.2xl0^^ cm~^ 
and Zi=A correspondingly. A 200 nm thick hydrogen-rich layer (ne=6xl0^^ cm~^) with 
transverse dimension of 2.5 fim is initially located at the back surface of the target. Two 
and three interaction stages have been simulated and the final proton energy (averaged 
over all protons) has been compared to that obtained in a single interaction scheme. The 
schematic diagram of multi-stage interaction setup is shown in figure [TJ In the multiple 
interaction scheme, the laser pulse of intensity Iq is split into n sub-pulses of equal intensity 
Iq/u that are made to interact with n targets. The proton layer is located only at the back 
surface of the first target. The remaining targets are devoid of any contaminant hydrogen- 
rich materials. 

In the two-stage setup, the proton layer gets accelerated by the electrostatic field devel- 
oped through the interaction of the first laser sub-pulse with the first substrate. The second 
laser sub-pulse in the mean time, travels to the second target, interacts with it and sets up 
a longitudinal electric field. The travelling proton layer passes through the second substrate 
and gets an extra boost from this electric field. If the arrival time for the second laser sub- 
pulse at the second target is properly adjusted, the proton layer gets an appreciable energy 
increase. In our numerical simulations we implement the time delay by performing several 
pre-simulation runs from which we determine the optimum timing of the laser pulse on every 
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subsequent target. It is interesting to note that the dependence of the final proton energy 
(and also the laser-to-proton energy transfer efficiency) on time delays exhibits a resonant 
behavior, so that any deviation from the most optimal distribution will lead to lower final 
proton energy. Therefore, the delay times will have to be properly adjusted in an actual ex- 
periment until the highest proton energy is achieved. Moreover, each individual interaction 
stage will introduce additional transverse divergence in the proton beam, which may lead to 
particle loss. The beam divergence is due to the generated transverse electric field (induced 
by the laser ponderomotive force gradient in the transverse direction). In order to limit this 
effect, one needs to assure that the laser pulse width is larger than the transverse size of the 
proton layer (e.g. by using layered— ctu..e. a. done. refe.e„ceH). 




The results of PIC simulations show that with the two-stage splitting the final average 
energy of the accelerated protons reaches Ep^^ =81.5 MeV, as opposed to Ep^^ =60.5 MeV 
(where the superscript denotes the number of interaction stages) for the conventional single 
target assembly, which is an increase of ~ 35%. Using the procedure described above, we 
have also simulated a 3-stage interaction scheme, in which case the main laser pulse is split 
into three sub-pulses of equal intensity I = Iq/3 that are made to interact with three targets 
with the same physical parameters described above. The final average proton energy in this 
3-stage setting reaches Ep = 96.5 MeV, which is ~ 60% energy increase as compared to 
the single interaction case or ~ 19% as compared to the 2-stage procedure. Figure [2] also 
shows the proton energy distributions for the three interaction stages. Gradual increase in 
the peak proton energy is clearly seen. 

As the number of splits n increases one should expect a gradual increase in the final proton 
energy. Of course there is a limit on the number of interaction stages, which would yield 
higher proton energies. This limit is related to the fact that the intensity of laser sub-pulses 
has to be high enough so that the laser ponderomotive force can still push electrons out of 
the target, thus setting up an accelerating electric field for protons. For estimation purposes, 
the number of splits n ~ Oq, where oq = eE/{'mcuj) is the laser relativistic parameter. In 
order to study how the proton energy changes with greater number of splits (n > 3) as well 
as splitting ratio between laser sub-pulses, we have used earlier developed model |m| for the 
longitudinal electric field, briefiy described in the next paragraph. 

The model is based on the fact that the accelerating electric field can be approximated by 
that of a charged cylinder of radius a and thickness 2ro, having the following mathematical 
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form (on the cylinder's 
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ro/|x| , Ixl > ro 



(1) 



where Qo is the charge of the target if all electrons are expelled, k = l/(47reo) and ri{t) is 



led electrons as a function of time that can be approximated by 



the proportion of the expe 
the following expression |l4j] 

f e-"(*-*«)' , t<to 

Vit) = 7 < , 2 

[ 5 + (l-5)e-'3(*-*o) ,t>to 

where 7 is the fraction of the electrons expelled at the peak of the laser pulse, S is the 
fraction of the initially expelled electrons that never return to the target, to is the arrival 
time of the peak of the laser pulse at the target, a = 41n2/r^ is a constant that depends on 
the pulse width (FWHM) used in the PIC simulation, (3 is the rate of return of the expelled 
electrons. These numerical factors are functions of laser intensity and have been tabulated 
using the PIC simulations. The equation of motion for a proton interacting with the field 
distribution ([1]) is: 

^x(t)=t;(t), (3) 

where rup is the proton mass, and e is the elementary charge. Eqs. ([3]) have been solved 
numerically for a wide range of splitting ratios x = h/Io {h is the intensity of the first laser 
sub-pulse) and a = h/H^ — x)h) {h is the intensity of the second laser sub-pulse) in the 
tree-stage interaction scheme. Figure [3] shows the final average proton energy as a function 
of the splitting ratio parameters normalized to the proton energy obtained from a single 
interaction stage. As one can see, the maximum in the proton energy occurs when x = 1/3 
and (7=1/2 (corresponding to three laser sub-pulses with equal intensities /o/3). Any other 
combination of splitting ratios leads to lower final proton energy. It should be noted that the 
two-stage results are recovered from this figure when one of the splitting parameters {x, o") 
are equal to or 1. In this case, the maximum in proton energy is reached at exactly equal 



splitting of the laser pulse into two sub-pulses with intensity /o/2. The fact that equal pulse 
splitting leads to maximum proton energy is not a mere coincidence and its origin will be 
discussed in the next section. Using Eqs. El we have also investigated how the final proton 
energy depends on the number of amplification stages, shown in figure HI As one can see, 
the final proton energy increases as the number of spitting stages becomes larger. 



III. INTERPRETATION OF OBTAINED RESULTS 

Results obtained in the previous section pose two interesting questions. Firstly, why 
splitting of the single interaction into multiple sub-stages gradually increases the amount 
of energy transferred from the laser pulse to protons? And secondly, why only equal laser 
pulse splitting leads to maximum proton energy? One possible answer to the first question 
can be obtained from the analysis of the microscopic dynamics of protons in the electric 
fields created through the interaction of each laser sub-pulse with each target. The analysis 
shows that the splitting procedure allows for more efficient utilization of the energy of the 
longitudinal electric fields in case of multiple interaction stages as opposed to that in the 
single, which gradually leads to higher proton energy. However, the same analysis does 
not yield a satisfactory answer to the second question. As mentioned in the introduction, 
proton acceleration by high power lasers can be qualitatively viewed from a thermodynamic 
perspective as an energy exchange process between hot (with initial temperature T/j) and 
cold (with initial temperature Tc) objects. Even though on a microscopic level the interaction 
potentials and the equations of motion for atoms on one hand and electrons and ions on the 
other are different, on a macroscopic level these two problems are quite similar in that they 
deal with energy flow and can be described using thermodynamic approach. The problem 
at hand may be formulated as follows: what is the most efficient method of exchanging the 
energy between the hot and cold objects, so that in the end the initially hot object becomes 
cold and initially cold object becomes hot? The most simple way by which one can exchange 
the heat between both objects is through putting them into the thermal contact with each 
other, so that in the end their final temperature is just half of the sum of their initial 
temperatures. The entropy change for this particular process corresponds to a maximum in 
the entropy gain, making it completely irreversible and least efficient in the sense of energy 



exchange between both objects. From a thermodynamic point of view 151], the efficiency of 
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the energy transfer from the hot object to the cold is at a maximum for those processes for 
which the entropy change tends to zero. Therefore, the problem is reduced to finding those 
processes that minimize the entropy gain. 

At this point, we would like to refer the reader to the Appendix section where we in- 
troduce the so-called "perfect heat exchange" problem 16|] and show that a very efficient 
heat exchange process between the hot and cold objects may be established using a splitting 
procedure. It is shown that just like in the present problem of the energy exchange between 
the laser pulse and protons, the splitting procedure in the "perfect heat exchange" problem 
also leads to more efficient energy exchange between the hot and cold objects. It is also con- 
firmed that the increased energy transfer efficiency is accompanied by decreasing increment 
in the total entropy change. This leads to the formulation of the following conjecture: the 
increased energy transfer efficiency between the laser pulse and protons in the multistage 
model is due to the reduction in the total entropy increment for the laser-target system. 

In addition, one can also arrive at the above conjecture with the help of PIC simulations. 
The total entropy change for the laser/target system can be estimated using the following 
expression, 

AS = ASi + AS, = AE (1/Te - l/T,) , (4) 

where AE is the energy transferred from the laser pulse to the target and is the initial 
temperature of plasma target (usually in the several keV range). Since the effective laser 
temperature Ti ^ (see below), one can neglect the term 1/T; in the brackets and for the 
multiple splitting case one obtains, 

1 " 

ASn = -y^AEjT,, (5) 

i=l 

where AEn is the energy transferred from the laser sub-pulse of intensity lo/n to the n*'' 
target. Using the PIC simulation results, we have calculated the total energy transferred 
from the laser pulse to the target (electrons, ions and protons), as a function of the laser 
splitting ratio I/Iq shown in Figure [5]a. As one can see, as the laser intensity becomes 
smaller, the amount of energy AEn transferred to the target is reduced, which eventually 
leads to a smaller increment in the entropy change for the total system. It is an expected 
result since the entropy production rate depends on the induced (by the laser pulse) tem- 
perature, density, pressure gradients in a plasma slab. Lower laser pulse intensity induces 
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smaller hydrodynamic gradients and lower heat fluxes, leading to smaller entropy production 
rates and correspondingly smaller total (particles and flelds) entropy change. At the same 
time, Figure Ob shows the efliciency of energy transfer from the laser pulse to protons as 
a function of laser energy obtained from the PIC simulations. As one can see, the energy 
transfer efficiency increases with decreasing amount of energy in the incoming laser pulse, 
confirming the earlier made conjecture concerning the relation between the energy transfer 
efficiency and the entropy change for the laser-target system. 

Finally, we would like to discuss the applicability of the thermodynamic interpretation 
of the multistage laser-proton acceleration. Expressions for the final temperatures or stored 
energies presented above have been derived under the assumption that both systems reach 
a state of thermodynamic equilibrium after the thermal contact has been established. In 
a case of laser interacting with the target, the thermal equilibrium between the laser and 
target is never achieved. Therefore, expression (]A3P cannot be directly used to evaluate 
the amount of energy stored in the proton beam. However, the fact that the thermal 
equilibrium between the laser and target is never achieved only results in lower value for 
the total amount of energy transferred between both systems, without changing the main 
conclusion that the efficiency of energy transfer from the laser pulse to protons increases 
with increasing number of splits. In addition, unlike in the "perfect heat exchange" problem 
where both the hot and cold reservoirs have been in their own corresponding equilibrium 
states (with well defined temperatures, entropies, etc.) prior to the thermal contact (the 
total closed system incorporating the hot and cold reservoirs however is a non-equilibrium 
system), in the laser-target system the laser pulse before the interaction with the target is 
not in a state of thermodynamic equilibrium. This fact however should not affect the main 
conjecture of the paper for the reasons presented below. 

In order to make a connection between the results of statistical mechanics as applied to 
the problem of heat exchange between hot and cold reservoirs on one hand and laser-target 
interaction on the other, definitions for the thermodynamic quantities such as temperature 
and entropy for a laser system are needed. One may also assume that the target substrate is 
in the state of thermodynamic equilibrium before and shortly after the interaction with the 
laser pulse, since the energy relaxation time r^e (or Maxwellization) between the hot and cold 



electron populations of the substrate is on the order of ~ 10 s 



231 ] for electron temperature 



of 1 keV, which is much shorter than the proton acceleration time (even though the relaxation 



10 



time for the ion component tu = Tee y ^ and that between electrons and ions r^i = T^e^ are 
longer than Tee, one may assume Maxwellian distribution for ions with its own temperature 
unequal to that for the electrons). The proton layer distribution function throughout the 
acceleration process may be assumed shifted Maxwellian with average velocity equal to 
that acquired by protons in the accelerating electric field (once the particles left effective 
acceleration region). Since the incoming laser light is not a black body radiation and is 
not in a thermodynamic equilibrium (in a sense that its entropy is not at a maximum), its 
effective temperature has certain ambiguity in its definition and still remains a matter of 



debate 2^ ■ Nonetheless, the thermodynamic definition of temperature given by expression, 

(6) 



1 - 

where S is the entropy of the laser pulse and E is its total energy, leads to the so called 



brightness temperature of the laser light 



calculated by counting states for identical bosons 



he entropy of the laser light is usually 
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24| . In this case the entropy is due to 



the randomness of the phase of each field quantum. Using expression ([6]) one can find the 



brightness temperature for laser light, given by the following 



15 



24| 



= MiTIM' ^ ^ 

where n is the mean photon occupation number, which is related to the laser specific intensity 
Icu (laser intensity per unit frequency) through the relation n = Ati^c^I^/Tiuj^ . Even for low 
intensity laser light, the mean occupation number is very large, so that the limit n » 1 
can be easily applied to the ultra-high intensity lasers. In this case, a very simple relation 
between the laser temperature (for mode uj) and its spectral intensity is obtained, 

ksT^ = rtX^I^. (8) 

This shows that each individual mode u has its own temperature T^. Note that for laser 
intensity Jq, carrying frequency ujq and spectral width Auj used in the PIC simulations, the 
laser temperature (at the central frequency c^o) reaches the staggering value of ~ 10^^ K. 
Thus, for any practical applications one can assume that the ultra-high intensity laser light 
is a system with infinite temperature. In this respect, the laser light can be considered to be 
a collection ot non-inte.acting (because of lineanty ot electrod^namic equations) =tatio,wy 
(see below) sub-systems, each with its own temperature [15^. Since different frequency 
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modes in the pulse do not interact with each other, the relaxation time for the propagating 
laser light into the equilibrium state is effectively infinite in the absence of matter. This in 
turn implies that a propagating laser pulse in itself may effectively be considered a static 
system in a state that resembles an equilibrium distribution (even though initially it was cre- 
ated by the system far away from the equilibrium and its spectral and angular distributions 
do not resemble the " black-body" radiation spectra) since its distribution function does not 
change in time as it propagates to and from the target. In addition, one may assume (ne- 
glecting the effects of the change in the laser angular and spectral distributions during the 
course of interaction) that the entropy change in the laser pulse as a result of its interaction 
with the target is smaller (in its absolute value) than that due to the target system (since 
laser effective temperature is infinite), so that the total entropy change is solely due to the 
target. As was mentioned in the introduction section, shorter relaxation times in the system 
lead to lower entropy increase and correspondingly higher energy transfer efficiency. Since 
the relaxation time for the electron component is equal the inverse of the electron-electron 
collision frequency Tee — ^eei'^Th)-! higher laser intensity leads to higher final electron tem- 
perature and correspondingly longer relaxation times and lower energy transfer efficiency 
(this argument also confirms the main conjecture of this work albeit from different consid- 
eration). In this respect, the arguments presented above give qualitative justification to the 
thermodynamic analogy between the " perfect heat exchange" problem on one hand and the 
multi-stage laser-plasma interaction model on the other. 

The main conclusion that one can draw from the example above is that splitting of the 
single interaction stage into multiple sub-stages is an effective way of reducing an irreversible 
component in the total interaction cycle no matter how this interaction looks like (laser- 
matter or matter-matter), thus increasing the effectiveness of the "pump". This is why the 
splitting procedure should also lead to higher proton energies in the laser-matter interaction 
experiments, since it increases the effectiveness of the energy transfer from the laser pulse 
to protons. 

IV. DISCUSSIONS AND CONCLUSION 

The problem of proton acceleration by high-power lasers has been revisited in this work 
using 2D PIC simulations as well as an analytical 3D model. It was shown that significant 
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energy gain in the final proton energy is possible if one introduces a multistage interaction 
scheme as opposed to a conventional single laser/target interaction setup. Many recent 
investigations concerning the proton acceleration looked at the kinematic/dynamic aspect of 
this problem, specifically the underlying physics behind the particle acceleration. However, 
the problem of how to increase the energy transfer efficiency from the laser to accelerated 
particles has not been addressed to the same degree of scrupulousness. As shown in the 
present work, the multistage interaction model offers significant gain in the efficiency of 
energy transfer from the laser to accelerated particles. A thermodynamic analogy has been 
offered to elucidate this effect. According to the model, the sphtting of a single interaction 
site into multiple stages is an effective way of reducing an irreversible component in the 
energy exchange process between the laser and protons. As a result more laser energy is 
transformed into proton kinetic energy. It was shown that in a three-stage setting, there 
is 60% increase in the energy efficiency of the laser accelerator as compared to a single 
interaction scheme. At the same time according to the results of our 3D model, it should 
be possible to increase the energy efficiency by more than 100% for a six-stage interaction 
setting without using more powerful lasers. Based on these results we conclude that the 
multi-staging procedure represents a step forward in increasing the energy efficiency of laser- 
ion accelerators with the potential of achieving significant increase in the final proton energies 
suitable for practical applications. 

APPENDIX A: PERFECT HEAT EXCHANGE PROBLEM AND ITS RELA- 
TION TO LASER-MATTER INTERACTION 

In what follows, we describe the "perfect heat exchange" problem and discuss its relation 
to laser-to-proton energy transfer efficiency. Let us suppose that we have hot and cold 
reservoirs (for a sake of simplicity we shall assume that both objects have the same size 
and mass and consist of an ideal gas) with initial temperatures and Tc correspondingly. 
The task is to find such a heat exchange process that would maximize the energy exchange 
between both reservoirs, so that in the end the hot reservoir becomes cold and the cold 
reservoir becomes hot. 

If initially hot and cold reservoirs are split into n equal pieces each and subsequently every 
individual hot piece is put into thermal contact with each individual cold piece (without 
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mixing them) in a sequential manner as shown in figure [6l then using the thermal balance 
equations, it can be easily shown that the final temperature of initially hot/cold objects 
(formed by putting back together initially hot/cold pieces to form new cold/hot objects) 
can be smaller/greater than [Th + T^) /2. The general expression for the final temperature of 
initially hot / cold reservoirs for n equal splittings is given by the following expressions [l^ : 

T, - - , (Al) 



rpifin) _ X]i=l i^i,nTh + 5n,iTc 



n 
j 



1 <^.ki^+J-k-l) 



-1)!^ 



2.+.(,_l)!Z-. (j-k)\ 



It can be shown that in the limit n ^ oo, T^-^*"^ = Tc and T^^^^ = Th and a perfect heat 
exchange process between hot and cold objects is established. Assuming that both objects 
are an ideal gas, the entropy change for the process involving n equal splittings has the 
following form: 

AS* nr=l i.^i,nTh + 5n,iTc) {5i,nTc + 5n,iTh) ^^^^ 



Cp 



nr=l {^i,nTh + 5n,iTc) {Si,nTc + S n,iTh 



where Cp is the specific heat capacity of the material. As the number of splits n increases, 
the entropy increment AS* decreases (remaining positive) and in the limit ^ oo, the 
entropy change AS* — 0. We would like to note that even though we used an ideal gas in 
the calculation of the entropy change, the same conclusion can be drawn if one were to use 
any other system. 

At this point we would like to draw the reader's attention to the fact that the multistage 
laser-target interaction scheme may be qualitatively viewed as a particular case of the "per- 
fect heat exchange" problem as presented above. In the proposed multistage model, a given 
laser sub-pulse interacts with a given substrate only once. However, the proton layer gets n 
energy kicks originating from the electric fields created by the interaction of the i-th laser 
sub- pulse with the ith. substrate for i G (1, n). In this respect the proposed multistage model 
may be viewed as a "heat exchange" problem between a small cold object of mass m and 
initial temperature (protons) and a large hot reservoir of mass M and initial temperature 
T/j (laser pulse), split into n small objects (sub-pulses) each having mass M/n. 

Moreover, in the derivation of eqs. lAll it is assumed that the interface between the hot 
and cold reservoirs is an ideal thermal conductor, so that there is no heat loss (e.g. radiative 
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losses) during the energy flow. In the laser-target interaction scheme the role of the interface 
is taken by the substrate material (electrons and ions of the substrate), which of course is 
not an ideal thermal conductor. This results in reduced energy gain for protons, while the 
splitting scheme remains beneficial. In addition, when the intensity of the laser pulse drops 
below a certain threshold value (so that electrons of the medium oscillate in phase with the 
incident light), the energy absorbed by the electrons of the substrate (through the absorption 
of the laser photons) gets reemitted back into the vacuum, so that in the end there is almost 
no net laser energy deposited in the target (classical skin effect). The threshold value for 
laser intensity at which this occurs is ao ~ 1. When the laser sub-pulse's amplitude drops 
below this value, there is virtually no net energy flow from the laser to the substrate and 
subsequently to the protons. Therefore, unlike the perfect heat exchange problem where 
the energy flow between the hot and cold reservoirs is always present no matter what the 
splitting ratio n is, in a case of laser-target interactions it would not be beneficial to split the 
laser pulse into too many sub-pulses. The threshold for the laser amplitude defined earlier 
determines the actual number of splitting stages for which one should expect energy fiow 
from the laser to protons. 

Using the thermal balance equations one can readily arrive at the expression for the final 
amount of energy stored in initially cold (with internal energy Ec) and hot (with internal 
energy Eh) objects as well as the corresponding entropy change, 

^F"^^ = (T^y^ + - ttAvt 1 , (A3) 



AS E, 



{fin) 



= Sin 



c 




Cp Ec 

where 5 = m/M. For a given value of 5 << 1, the final energy stored in a small object 
E^^"'^ increases with the number of splits n, quickly approaching the limiting value E^S. 
The total entropy change also decreases with the number of splits as in the general case 
discussed earlier, however unlike in the "perfect heat exchange" problem, it never goes to 
zero (in the limit n — > oo) but approaches a finite value, which depends on 5. This indicates 
that there is a limit on the efficiency of the energy transfer between the two objects in this 
particular case, leading us to conclusion that there may also be a limit on how much of the 
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laser energy can be transferred to protons in the multistage laser-matter interaction scheme. 
This limit can also be seen in Figure H] where the final proton energy gradually approaches 
saturation with the number of interaction stages. The same behavior was also observed in 
the PIC simulation results. In addition, the PIC simulation also showed that only equal 
laser pulse splitting led to the highest proton energy. Any other combination of splitting 
parameters would lead to lower proton energy. As we mentioned earlier, this peculiarity is not 
coincidental. Analyzing the thermal balance equations in the case of arbitrary distribution 
of the masses of hot small objects (energies of laser sub- pulses), it can be shown using the 
variational principle that only homogeneous distribution (equal splitting) will lead to the 
maximum in the final temperature (energy of protons) of initially cold object. 
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FIG. 1: Proposed experimental setup for increased proton energy; a) - conventional double-layer 
target geometry; b) - two-stage proton acceleration; c)-three-stage proton acceleration 
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FIG. 2: Proton energy distributions for three different interaction stages. Solid line represents a 
single interaction stage (no laser splitting), dotted line represents a double interaction stage (single 
laser splitting) and dashed line represents a tripple interaction stage (double laser splitting) . 
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FIG. 4: Final proton energy as a function of the number of amplification stages 



21 




22 



FIG. 5: a) Energy transferred from the laser sub-pulse to the target versus the laser splitting ratio 
n = Iq/I. b) Normalized energy transfer efficiency from laser pulse to protons versus normalized 
laser intensity, /q = 1.9 x 10^^ W/cm^ 
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FIG. 6: Schematic diagram for two heat exchange processes, a) The hot and cold reservoirs are 
put into thermal contact with each other leading to temperature equalization. The entropy gain is 
maximal for this process, b) The hot and cold reservoirs are split into n pieces each that are put 

into thermal contact with each other in a sequential manner. The limit n oo, corresponds to 
totaly reversible heat eachange process with zero entropy gain. 
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